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Reports  and  Publications 

Progress  Reports,  periodically  through  duration  of  research  effort 
Final  Report,  this  document 

Research  Paper,  "Line  Mixing  in  the  v  -Band  of  C0„,"  to  be  submitted 

2  ^ 

for  publication  in  Applied  Optics 


Objectives 

The  original  objectives  of  this  research  project  were: 

1)  to  obtain  experimental,  high-resolution  absorption  spectra 

of  C0„  in  the  vicinity  of  the  v  -band  near  667  cm-1. 

2  2 

2)  to  include  line  mixing  effects  into  the  theoretical  description 
of  the  absorption  spectrum. 

In  order  to  meet  the  first  objective  a  tunable  diode  laser  spectrometer 
was  acquired  from  Laser  Analytics,  Inc  -  spectrometer  was  installed, 

aligned  and  calibrated.  The  necessary  vacuum  and  gas-handling  apparatus 
was  constructed  and  tested  to  allow  absorption  spectra  of  research  grade 
gaseous  CO^  to  be  taken  at  known  pressures.  In  order  to  analyze  the 
absorption  spectra  the  experiment  was  placed  on  line  to  a  PDP  11/60  computer 


located  in  an  adjacent  room  in  the  physics  building  at  New  Mexico  State 

ii 


University.  The  necessary  computer  software  was  written  to  permit  analysis 

of  experimental  absorption  spectra.  Preliminary  absorption  spectra  were 

obtained  for  CO.,  at  low  pressures  ~  1  Torr  (the  Doppler-limited  regime) . 

Similar  data  had  previously  been  reported  by  Aronson  et  alJ  and  has  since 

8  9 

been  reported  in  a  series  of  papers  by  Planet  and  Tettemer  ’  (the  reference 

numbers  refer  to  the  bibliography  at  the  end  of  this  report) .  The  preliminary 

absorption  spectra  at  low  pressures  were  intended  to  serve  as  a  calibration 

of  the  spectrometer  resolution  and  sensitivity.  The  intention  was  to  follow 

this  preliminary  series  of  low  pressure  experiments  with  experiments  at 

higher  pressures  where  line  mixing  effects  could  be  studied.  The  results 

of  the  work  described  in  this  paragraph  were  reported  to  ARO  in  4  Progress 

Reports  covering  the  period  April,  1977  -  April,  1979. 

At  this  point,  as  noted  in  a  Progress  Report  sent  to  ARO  in  the  fall 

of  1979,  the  laser  diode  failed.  A  replacement  diode  was  ordered  but 

difficulties  in  the  fabrication  of  diodes  at  the  desired  frequency  of 

667  cm-1  resulted  in  an  excessive  delay  before  delivery  of  the  diode.  It 

was  realized  that  alternative  sources  of  diodes  should  be  sought.  One 

source  was  found,  the  Atmospheric  Sciences  Laboratory  (ASL)  at  White  Sands 

Missile  Range,  ASL  kindly  loaned  us  3  diodes  that  covered  the  range  from 

950-1150  cm-1.  Since  investigating  the  v  -band  of  CO,  near  667  cm  1  was 

2  z 

impossible  with  these  diodes  alternative  spectral  problems  were  considered. 

An  initial  group  of  candidates  were  the  CO^  laser  bands  at  950  and  1050  cm-1 
but  both  preliminary  absorption  measurements  and  absorption  calculations 
indicated  that  absorption  from  these  bands  was  too  weak  to  be  detectable 
with  our  current  spectrometer.  The  initial  designs  were  being  drawn  for 


iii 


a  longer  path  length  absorption  cell,  which  would  increase  the  sensitivity 
of  our  spectrometer,  when  the  cryogenic  compressor  in  the  closed  cycle 
helium  refrigerator  developed  a  vacuum  leak.  At  this  time  the  post-doctoral 
fellow,  who  had  been  working  on  the  project  with  me,  was  forced  by  personal 
reasons  to  leave  the  project.  I  was  able  to  hire  a  graduate  student 
half-time  for  an  entire  academic  year  by  using  the  remainder  of  the  post¬ 
doctoral  fellowT's  salary.  The  results  of  the  efforts  outlined  in  this 
paragraph  we  described  in  Progress  Reports  to  ARO  in  the  spring  and  fall  of 
1980,  and  the  personnel  action  which  resulted  in  the  hiring  of  the  graduate 
student  was  outlined  in  a  letter  to  ARO  in  September  of  1980. 

Through  the  dilligent  efforts  of  the  graduate  assistant  the  cryogenic 
compressor  was  placed  back  in  working  condition  in  the  early  spring  of  1981. 

At  this  time  extensive  tests  revealed  the  unfortunate  fact  that  all  diodes 
in  our  possession,  the  three  diodes  on  loan  from  ASL  and  the  replacement 
diode  purchased  over  a  year  before,  were  either  inoperative  or  had  insufficient 
power  levels  for  practical  spectroscopic  work.  In  particular  the  output 
power  level  of  the  667  cm-1  diode  was  much  too  low  to  permit  the  higher 
pressure  experiments  to  be  done  on  the  v  -band  of  C09 .  This  unfortunate 
occurrence  thus  deprives  this  research  effort  from  experimental  data  on  the 
effect  of  line  mixing  in  this  theoretically  and  practically  significant  region 
of  the  spectrum. 

The  second  objective  of  this  research  work  has  been  achieved.  A 

theoretical  description  of  the  absorption  spectrum  of  the  v  -band  of  CO 
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has  been  obtained  which  includes  the  effect  of  line  mixing.  The  results 
of  this  work  are  given  in  the  following  report.  An  edited  version  of 
this  report  will  be  submitted  shortly  to  Applied  Optics  for  publication. 
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ABSTRACT 


A  theoretical  description  of  line  mixing  in  the  Q-branch  of  an 
infrared  rotation- vibration  band  is  obtained  in  the  weak  coupling 
approximation.  The  resulting  absorption  band  shape  is  applied  to  the 
case  of  the  v^-band  of  CO^  near  667  cm-1.  Theoretical  absorption  profiles 
are  illustrated  in  selected  spectral  regions  for  several  values  of  the 
pressure . 


I.  Introduction 


The  determination  of  atmospheric  absorption  is  important  in  the  analysis 

of  a  variety  of  current  scientific  and  technological  problems  as,  for 

example,  in  the  use  of  satellite-born  infrared  radiometers  to  retrieve 

1  2 

atmospheric  temperature  profiles.  In  the  infrared  region  of  the  spectrum 

the  shape  of  an  absorption  band  is  often  modelled  by  summing  the  absorption 

at  each  frequency  due  to  the  individual  lines  in  the  band.  Line  parameter 

data  is  conveniently  available  for  many  molecules  of  atmospheric  interest 

3 

from  the  AFGL  compilation.  The  absorption  profiles  near  line  center  are 
nearly  Lorentzian  at  higher  pressures:  at  lower  pressures  Doppler  broadening 
effects  are  also  significant. 

The  above  approach  is  adequate  at  pressures  where  there  is  no  appreciable 
overlap  between  lines  in  the  absorption  band.  As  the  pressure  is  raised 
line  mixing  effects  become  significant  and  cause  the  absorption  to  depart 
from  the  summation  over  individual  lines  described  above.  Line  mixing  effects 
have  been  observed  in  infrared^’'’  and  Raman  spectra*^  at  pressures  in  excess 
of  one  atmosphere.  Other  molecular  systems,  such  as  the  Q-branches  of  infrared 
and  Raman  bands  and  the  rotational  Raman  spectra  of  symmetric  top  or  slightly 
asymmetric  top  molecules,  have  groups  of  overlapping  lines  even  at  the  moderate 
pressures  characteristic  of  the  earth's  atmosphere:  line  mixing  should  then 
be  exhibited  by  these  systems  at  pressures  below  one  atmosphere. 

In  this  paper  we  consider  line  mixing  in  the  Q-branch  of  the  v^-band 

1  ° 

of  C0o  near  667  cm-1.  This  band  is  important  in  radiometric  applications, 

7  8  9 

has  been  studied  extensively  at  low  pressures  ’  ’  and  line  parameter  data 

3  7  8  9 

are  available  from  several  sources.  ’  ’  ’  The  Q-branch  consists  of  a  series 
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of  closely  spaced  lines,  split  by  the  rotation-vibration  interaction.  A 

qualitative  interpretation  of  the  pressure  dependence  of  the  Q-branch  band 

shape  may  be  given  following  the  discussion  of  Alekseyev  et  al.^  At  low 

pressure,  where  the  rotationallv  inelastic  collision  frequency  is  much  less 

than  the  rotation-vibration  interaction,  the  lines  are  well  resolved:  in 

this  regime  individual  lines  are  broadened  by  collisions.  As  the  pressure 

is  increased  the  lines  began  to  overlap  and  collisions  broaden  and  mix  the 

lines  creating  interference  terms  in  the  band  shape.  Finally,  at  high 

pressures,  where  the  rotationally  inelastic  collision  frequency  is  large 

compared  to  the  rotation-vibration  interaction,  the  lines  in  the  band  merge. 

In  this  regime  modifications  in  the  rotation-vibration  interaction  result  in 

all  lines  in  the  band  having  the  same  vibration  frequency:  the  band  shape 

narrows  with  increasing  pressure,  producing  the  phenomenon  of  motional 

narrowing.  For  the  present  case  of  the  Q-branch  of  CO^  we  may  obtain  a 

measure  of  the  rotation-vibration  interaction  by  considering  the  line  splitting 

3 

in  the  neighborhood  of  the  intensity  maximum.  These  lines  are  split  by 
about  .068  cm-1  and  it  may  then  be  seen  that,  for  T  -  300  K,  the  rotationally 
inelastic  collision  frequency  becomes  equal  to  the  rotation-vibration  inter¬ 
action  at  a  pressure  -690  torr.  Thus,  line  mixing  should  be  significant  in 
the  Q-branch  at  pressures  below  one  atmosphere  but  motional  narrowing  should 
not  be  significant  until  pressures  of  several  atmospheres  are  reached. 

The  determination  of  the  band  shape  in  the  presence  of  line  mixing 
presumes  a  knowledge  of  the  elements  of  the  transition  rate  matrix.^  Direct 
calculation  of  the  elements  of  this  matrix  requires  that,  for  an  assumed 
inter-molecular  potential,  an  average  be  taken  over  the  dynamical  variables 
describing  a  collision.  This  is  a  formidable  numerical  problem  which,  in 
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principle,  must  be  solved  for  each  value  of  temperature  appropriate  to  the 

application.  To  circumvent  this  difficulty  we  assume  in  the  remainder  of 

this  paper  that  the  significant  collisions  which  determine  the  band  shape 

12 

satisfy  the  weak  coupling  approximation.  In  this  approximation  the  rate 

matrix  is  dominated  by  its  diagonal  elements.  In  the  present  case  we 

consider  only  the  dominant  isotopic  species  i2C^09  for  which  symmetry 

considerations  require  J  to  be  even.  This  implies  that  in  most  collisions 

the  angular  momentum  either  remains  the  same  (AJ  =  0)  or  changes  at  most  bv 

two  units  (AJ  =  ±2) .  The  weak  coupling  approximation  lias  been  used  by 
13 

Rosenkranz  to  adequately  model  the  microwave  spectrum  of  atmospheric  oxygen 

4 

near  5  mm.  The  work  of  Gordon  and  McGinnis  on  the  infrared  spectrum  of  CO 
near  2150  cm-1  also  indicates  the  dominance  of  the  near-diagonal  terns  in 
the  rate  matrix. 


II.  Weak  Coupling  Approximation 

14 

In  the  Heisenberg  picture  of  spectroscopy  the  absorption  band  shape  is 
determined  by  the  Fourier  transform  of  the  transition  dipole  moment  correlation 
function : 


<h(0)*ij(t)>  =  d-exp[i(«o+  i  31 )  1  *  p  -  d  ,  (1) 

where  d  and  p  are  vectors  of  intrinsic  line  amplitudes  and  fractional  popu¬ 
lations,  respectively,  is  a  matrix  of  line  frequencies  and  II  is  the 
transition  rate  matrix.  The  brackets  indicate  a  thermal  average.  Diagonal 
elements  of  the  rate  matrix,  K  ,  give  the  rate  at  which  collisions  transfer 
amplitude  away  from  line  k  to  all  other  lines  in  the  spectrum  while  the 
off-diagonal  elements,  Rive  minus  the  rate  at  which  collisions  transfer 


amplitude  from  line  k  to  line  i.  The  rate  matrix  may  be  conveniently  expres; 


4 


in  terms  of  a  cross  section  matrix  o  as: 


f!  =  nva  , 


(2) 


where  n  is  the  absorber  number  density,  and  v  is  the  average  absorber- 
perturber  relative  speed.  The  cross  section  matrix  has  elements: 


Ik 


1 

v 


<v(5fk  '  Sik)>2"bdb' 


Ci) 


where  b  is  the  impact  parameter  and  S  is  an  element  of  the.  scatterin’--  t:.i *  i' i 

k 

that  gives  the  fractional  amplitude  transferred  by  a  collision  into  1 ine 
from  a  unit  amplitude  in  line  k  before  the  collision. 

Gordon"**  has  developed  a  semi-classical  descr  iption  of  me. ecu :  nr 
collisions  in  order  to  simplify  the  problem  of  the  determination  of  the 
scattering  matrix  elements ,  S^.  In  this  description  the  frequencies  •  f  inc¬ 
lines  in  the  spectrum  are  assigned  to  appropriate  classical  rotation  fre¬ 
quencies  by  means  of  the  correspondence  principle.  Thus,  fur  an  infrared 
rotation-vibration  band  the  P-  and  R-b ranch  lines  arise  from  rotation  of  the 
transition  dipole  moment,  in  opposite  senses,  at  frequencies  appropriate  to 
the  overall  rotation  of  the  molecule.  The  Q-br.inch  linos,  however,  arise 
from  quantum  mechanical  transitions  between  different  vibrational  annular 
momentum  states.  The  corresponding  classical  rotational  frequencies  f or 
the  Q-branch  lines  are  then  much  higher  than  those  of  the  P  and  ’{-branch 
lines.  This  case,  Hund's  case  (a),  has  been  considered  by  Gordon*'  who 
finds  that  in  the  classical  description  collisions  do  not  couple  the  ‘  branch 
to  either  the  P  or  the  R-branches  (the  P  and  R-b  ranches  are  collisior.ally 
coupled  principally  by  the  so-called  "reorientation  effect"  which  is  responsi: 
for  the  collapse  of  the  whole  band  to  a  single  lino  at  extremely  high  piaosurc 


but  this  coupling  does  not  concern  us  here).  Thus,  in  the  study  of  the 

Q-branch  absorption  we  need  consider  only  the  effect  of  collisions  that 

transfer  amplitude  between  different  0-branch  lines. 

13 

In  the  weak  coupling  approximation  the  assumption  is  made  that  in  most 

collisions  the  angular  momentum  either  remains  the  same  or  changes  at  most 

by  two  units.  Thus  each  column  of  the  rate  matrix  contains:  a)  the  diagonal 

element  IT^;  2)  an  element  coupling  line  k  to  the  adjacent  line  at 

lower  frequency;  and  3)  an  element  II  „  coupling  line  k  to  the  adjacent  line 

k*T*2  k 

at  higher  frequency.  Because  of  the  definitions  of  the  elements  of  the  rate 
matrix,  conservation  of  probability  requires  that  each  column  satisfy: 


1  n 


ik 


0  , 


(4) 


which  in  the  weak  coupling  approximation  becomes; 


k-2k 


\k  + 


“k+2k  ° 


(■r» 


There  is  a  second,  rigorous,  relationship  between  elements  of  the  rate  matrix 
required  by  the  principle  of  detailed  balance.^  This  requirement,  which 
insures  that  the  absorbing  system  remains  in  thermal  equilibrium,  results  in 
elements  of  the  rate  matrix  satisfying  the  relation: 


Pk  \-2k  “  Pk-2 


'kk- 


(6) 


The  diagonal  elements  of  the  rate  matrix  are  simply  the  isolated  line- 

3  8  9 

widths  and  are  available  from  several  sources.  ’  ’  Equations  (5)  and  (6) 
therefore  suffice  to  determine  all  elements  of  the  rate  matrix  in  the  weak 


coupling  approximation. 


6 


Following  Rosenkranz  we.  solve  for  the  elements  in  descending  order. 

Hie  reason  for  this  is  that  Eq.  (5)  is  only  an  approximate  relation,  obtained 
by  truncating  Eq.  (4),  whereas  Eq.  (6)  is  exact.  Thus,  when  we  apply 
Eq.  (5)  to  the  highest  frequency  line  conside.red  we  set  'T  _  =  0  and  introduce 
a  large  relative  error  into  the  values  of  the.  elements  for  that  line.  But, 
this  error,  and  hence  the  error  introduced  into  the  rate  matrix  elements,  is 
damped  by  successive  applications  of  Eq.  (6).  Provided  we  select  a  maximum 
line  frequency  for  which  the  Boltzmann  factor  is  extremely  small,  the  initial 
error  will  produce  only  a  small  effect  in  lines  for  which  the  population 
is  large. 

In  the  present  case  we  have  included  lines  up  to  .1  =  106;  data  for 

3 

these  lines  are  available  from  the  AFGL  compilation.  We  thus  set: 

T  =  0  .  (7) 

108,106 

Application  of  Eqs.  (5)  and  (6)  then  gives: 

=  r  (  P106-2  (m-i) 

106-2 (nri-1) ,106- 2m  A  '  106-2 (m-i)  p.„,  0 

(B) 

where  the  linewidth  a  =  R  is  the  diagonal  element  for  the  m^1  column  of 

m  mm 

the  rate  matrix.  Equation  (5)  may  be  used  to  generate  the  corresponding  rate 

matrix  elements  coupling  each  line  to  the  adjacent  line  at  higher  frequency, 

except  for  the  lowest  frequency,  J  =  2,  line  which  must  be  handled  separately. 

The  reason  for  this  is  again  the  result  of  the  approximate  nature  of  F.q .  (5). 

For  the  J  =  2  line  we  must  set  n  =  0.  In  this  case  Eqs.  (5)  and  (6)  may 

02 

not  be  satisfied  simultaneously.  We  choose  to  satisfy  the  rigorous  condition, 
Eq.  (6),  which  in  this  case  takes  the  form: 
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>  fc 


p  h  =  p  n  .  (9) 

2  42  4  24 

The  approximate  relation,  Eq.  (5),  is  then  not  satisfied  for  the  lowest 
frequency  line. 

The  absorption  coefficient  may  be  computed  from  the  Fourier  transform  of 
the  transition  dipole  correlation  function,  Eq .  (1).  The  result  is : 


Y,  x  _  8m3nw(l-e  hcw'/kT) 

K(m) - 3hc -  F(U)  > 


(10) 


where 


F(u)  =  —  Im  d-(wl  -  to  -  ill)-1*p*d  , 
t  0 


(11) 


with  1  the  unit  matrix.  The  formal  mathematical  problem  in  obtaining  a  solu¬ 
tion  to  Eq.  (11)  is  to  diagonalize  the  non-Hermitian  matrix  (u0  +  ill)  .  A 
generalized  approach,  the  QR  algorithm,^  exists  for  performing  this  diagonal- 

ization,  as  well  as  an  approximate  perturbation  theory  solution.  Following 
13 

Rosenkranz  we  elect  to  use  the  perturbation  approach.  It  is  shown  in 
Appendix  A  that,  to  first-order  in  the  elements  11^  (which  is  equivalent  to  a 
first-order  expansion  in  the  pressure),  the  eigenvalues  of  (wq  +  ill)  are: 


Xk  =  “k  +  iHk  * 


(12) 


while  the  eigenvectors  have  elements: 


Xik  =  1 


n 


i  =  k  (13a) 

«  i  ~ — ,  i  /  k  (13b) 

w,  -w„ 
k  l 

Then  it  may  be  shown,  using  Eqs .  (6,12,13),  that  the  following  expression 


resul ts  for  F (w) : 


In  the  following  section  we  will  apply  these  results  to  the  analysis  of  the 
Q-branch  band  shape. 


Two  approximations  were  made  in  the  derivations  of  Eqs .  (14)  and  (15); 
the  weak-coupling  approximation  and  the  linear  approximation  to  the  eigen¬ 
values  and  eigenvectors  of  the  matrix  (w  +  ill)  .  In  order  to  disentangle  the 

0 

effects  of  these  two  approximations  on  the  computed  band  shape  in  a  simple 

case,  we  have  considered  a  synthetic  spectrum  containing  just  4  lines.  Assumin 

the  weak  coupling  approximation  to  be  valid  we  have,  with  the  aid  of  Fqs .  (5) 

and  (6),  rigorously  diagonalized  the  matrix  (w  +  ill).  The  details  of  this 

0 

calculation  are  given  in  Appendix  B.  Although  the  synthetic  4-line  spectrum 
is  not  equivalent  to  the  0-branch  spectrum  at  elevated  temperatures  (e.g ., 
at  T  -  300  K  the  maximum  Q-branch  intensity  occurs  for  J  =  16,  the  line 
in  the  band) ,  at  sufficiently  low  temperatures  and  pressures  the  synthetic 
spectrum  will  approximate  the  Q-branch  spectrum.  In  the  following,  section  we 
investigate,  for  low  temperatures  and  pressures,  the  correspondence  between 
the  synthetic  4-line  spectrum  and  the  spectrum  obtained  using  eqs.  (]4)  and 
(15),  and  hence  get  a  measure  of  the  validity  of  the  linear  approximation. 


J 
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III.  Absorption  Band  Shape 

We  consider  the  absorption  coefficient  K(w)  ,  given  by 
the  spectral  function  F(aj)  is  given  by  Eqs .  (14)  and  (15). 
dR  may  be  expressed  in  terms  of  rotational  line  amplitudes 


Eq .  (10),  where 

The  line  amplitude 
3 

D  according  to: 
k 


\  \  •  06) 

where  $v  is  a  reduced  vibrational  matrix  element.  It  is  convenient  to 
define  an  integrated  band  intensity  S°  by:^ 

32rJ  a>  (1-e  v  ) 

„  o  V  |  „  .  9 


where  w  is  the  frequence  of  the  band  origin.  Then  in  terras  of  S°  we  mav 
v  v 

express  the  absorption  coefficient,  neglecting  a  small  non-rigidity  factor, 


as  : 


o  .  -hcwA-T, 
nS  u>(l_e  ) 

K(w)  =  - - - ~ /pep -  G(w)  , 


u  (1-e  hcwv/kT 


) 


(18) 


where  the  spectral  function  G(u)  is  defined  by: 


G(u>)  =  Gp(u)  +  GR(U)  +  Gq(u).  0  9) 

In  Eq.  (19)  the  spectral  functions  Gp(u),  GR(u))  and  0^(w)  are  defined  in 

terms  of  the  rotational  line  line  amplitudes  D  .  For  the  case  of  the  v 

k 

1 8 

band  of  CO^  we  have: 
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a 

II 

f2 (k-i) y 

1  2k+l  J 

(20a) 

l .. 

<■ 

r 2  (k+2) , ; 

1  2k+l  J 

(20b) 

■2- 

[2]4'  . 

(20c) 

Then  we  have : 


Vw> 


I  V  A  f  k_1  1  V 

TT  £  Pkl2(2k+1)  J(w-M  F  +  uf  ’ 

k  k  k 


(21a) 


gr(“> 


_  1  Y  r  k+2  , _ ak 

_  tt  ^  Pkl2(2k+2)  J  (u-uk)z  +  a£  ’ 


(21b) 


where 


Vw> 


«k  4-  (u-a,k)Yk 
(u-U)k)  2  +  a£  ’ 


Y 


k 


Jtk 


u,  -u) 

k 


l 


(21  c) 


(22) 


Equation  (18)  for  the  absorption  coefficient,  together  with  the  definitions 
of  the  spectral  functions  given  by  Eqs .  (21a)-c)  and  (22),  may  then  be  used 
to  compute  an  absorption  band  shape.  The  P  and  R-branch  lines  which  overlap 
the  Q-branch  must  be  included  in  order  to  get  a  meaningful  band  shape,  although 
there  is  no  mixine  between  them.  Figures  1-6  are  plots  of  the  computed 
absorption  for  two  spectral  regions;  below  the  band  origin  from  661-667  cm  1 , 
and  the  region  that  includes  the  first  4  Q-branch  lines  from  667.366-667.469 
cm*1.  Pressures  of  10,  50,  and  100  torr  are  considered,  all  for  a  temperature 
of  300  K. 

The  plots  of  the  661-667  cm-1  region  illustrate  the  low  frequency  wing  of 
the  Q-branch  including  the  first  four  P-branch  lines.  The  solid  curve  is  the 


FIG.  1  Theoretical  absorption  coefficient  versus  frequency, 

661-667  cm-1  region,  v^-band  of  CO.,;  pressure  10  Torr; 

temperature  300  K;  solid  curve  with  line  mixing;  unconnect 
points  with  no  line  mixing. 


absorption 
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\  FIG.  2  Theoretical  absorption  coefficient  versus  frequency, 

!  661-667  cm'1  region,  v2~band  of  CO.,;  pressure  50  Torr; 

;  temperature  300  K;  solid  curve  with  line  mixing;  unconnct 

’  points  with  no  line,  mixing.. 


ABSORPTION 


Theoretical  absorption  coefficient  versus  frequency, 
661-667  cm-1  region,  v  -band  of  C0~;  pressure  100  Torr; 
temperature  300  K;  solid  curve  with  line  mixing;  unconnec 
points  with  no  line  mixing. 


ABSORPTION 


Theoretical  absorption  coefficient  versus  frequency, 
667.366-667.469  cm-1  region,  v2~band  of  CO^;  pressure 

10  Torr;  temperature  300  K;  solid  curve  with  line  mixing 
unconnected  points  with  no  line  mixing. 


FIG.  5  Theoretical  absorption  coefficient  versus  frequency, 
667.366-667.469  cm-1  region,  v^band  of  CO^;  pressure 

50  Torr;  temperature  300  K;  solid  curve  with  line  mixing; 
unconnected  points  with  no  line  mixing. 


FIC.  6  Theoretical  absorption  coefficient  versus  frequency, 
667.366-667.469  cm-1  region,  v^-band  of  C02;  pressure 

100  Torr;  temperature  300  K;  solid  curve  with  line  mixing; 
unconnected  points  with  no  line  mixing. 
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computed  absorption  with  line  mixing  included  while  the  unconnected  points 
illustrate  the  case  of  no  line  mixing  (all  Y  =  0  in  Ec .  (22)).  The  effect 
of  line  mixing  is  seen  to  be  mosL  apparent  in  the  troughs  between  the 
P-branch  lines,  reducing  the  absorption  there  by  greater  than  a  factor  of  10 
just  below  the  band  origin  at  667  cm”1.  The  effect  of  line  mixing  is  also 
observed  to  be  linearly  dependent  on  pressure  in  the  troughs,  appropriate  to 
the  first-order  approximation  made  in  inverting,  the  matrix  (u  +  if) .  Truncation 
of  the  line  is  caused  by  a  lack  of  resolution  in  the  plotting  program. 

The  plots  of  the  spectral  region  near  667  cm-1  illustrate  the  dependence 
of  the  line  mixing  effect  on  the  splitting  of  the  lines. ^  The  rotation- 
vibration  interaction  results  in  an  increased  line  splitting  with  increasing  -J, 
with  the  smallest  splitting,  .015  cm-1,  occurring  between  the  lowest  members 
of  the  Q-branch,  and  Q^.  The  effect  of  line  mixing  should  then  be  first 
evident  in  the  lowest  Q-branch  lines,  progressing  with  increasing  pressure  to 
higher  Q-branch  lines.  This  is  clearly  shown  in  Figs.  4-6.  At  the  lowest 
pressure  of  10  torr  line  mixing  has  already  reduced  the  absorption  a  factor  ~2 
below  Q2,  whereas  the  absorption  near  Qg  remains  virtually  unaffected.  At 
the  higher  pressures  of  50  and  100  torr  line  mixing  is  seen  to  have  a  progress¬ 
ively  greater  effect  on  the  higher  Q-branch  lines. 

Figure  7  illustrates  the  effect  of  the  linear  pressure  approximation 
made  in  obtaining  the  eigenvalues  and  eigenvectors  of  the  matrix  (uQ  +  ill)  . 

The  solid  line  gives  the  calculated  absorption  for  a  synthetic  4-line  Q-branch. 
The  weak-coupling  approximation  was  used  to  compute  this  spectrum  but  the 
matrix  (ojq  +  ill),  as  shown  in  Appendix  B,  was  rigorously  diagonalized.  The 
unconnected  points  represent  the  calculated  absorption  in  the  linear  pressure 
approximation.  The  synthetic  4-line  spectrum  does  not  closely  approximate 


Synthetic  absorption  coefficient  versus  frequency, 

667 . 366-667 .469  era-1  region,  v2~band  of  CO^;  pressure 

10  Torr,  temperature  25  K;  solid  curve,  exact  4-line 
spectrum;  unconnected  points,  linear  pressure  approximation. 
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the  Q-branch  at  elevated  temperatures  but  at  sufficiently  low  temperatures 
the  population  of  the  rotational  levels  falls  off  quickly  and  the  synthetic 
four-line  spectrum  becomes  an  increasingly  useful  approximation.  Figure  7 
illustrates  the  case  of  T  =  25  K,  for  which  the  peak  population  occurs  for 
the  line  (note  that  this  represents  a  "non-physical"  case  since  gaseous 
C0?  could  not  exist  at  this  temperature  in  thermal  equilibrium.  The  purpose 
here,  however,  is  simply  to  compare  two  theoretical  descriptions  of  the 
absorption,  and  the  non-physical  nature  of  the  temperature  is  of  no  significance). 
The  major  discrepancy  occurs  in  the  vicinity  of  the  line,  below  line  center, 
and  in  the  trough  between  the  Q2  and  0^  lines,  reaching  a  maximum  value  of 
about  20%  there.  In  the  remainder  of  the  spectral  region  the  two  computed 
spectra  are  virtually  indistinguishable.  A  comparison  of  Fig.  7  with  Figs.  1-6 
reveals  that  the  exact  calculated  absorption  enhances  the  effect  of  line 
mixing  as  compared  with  the  linear  approximation.  For  example,  below  the 
center  of  the  Q2  line,  the  exact  absorption  falls  below  the  linear  result, 
whereas  in  the  trough  between  the  Q2  and  lines  the  exact  absorption  results 
in  the  depth  of  the  trough  being  reduced  further  than  is  predicted  by  the 
linear  theory.  We  conclude  that  the  effect  of  line  mixing  will  be  at  least 
as  large  as  the  prediction  of  the  linear  pressure  theory. 

IV.  Conclusions  and  Recommendations 

In  this  report  we  have  computed  the  effect  of  line  mixing  on  the  Q-branch 
of  the  v  -band  of  C0„  near  667  cm  1 .  The  weak  coupling  approximation  was 

i  1 

used  to  generate  the  elements  of  the  transition  rate  matrix.  The  spectrum 
was  computed  by  means  of  a  linear  approximation,  which  was  used  to  diagonalize 
the  matrix  (u>  +  in) .  Absorption  band  shapes  were  computed  for  two  spectral 


23 


regions,  below  the  band  origin  from  661-667  cm-1  and  in  the  region  of  the 
first  4  Q-branch  lines  above  667  cm  J.  Finally,  the  absorption  spectrum  using 
the  linear  approximation  was  compared  with  a  synthetic  spectrum  in  which 

the  matrix  (w  +  ill)  was  rigorously  diagonalized. 

0 

In  the  linear  pressure  approximation  line  mixing  is  seen  from  Figs.  1-6 
to  have  a  significant  effect  on  the  computed  spectrum.  For  the  pressures 
considered,  between  10  and  100  torr,  the  effect  of  line  mixing  is  largest  in 
the  region  just  below  the  band  origin  near  667  cm_i  and  in  the  troughs  between 
the  lowest  Q-branch  lines,  around  667.4  cm-1.  We  may  conclude  from  the  results 
illustrated  in  Fig.  7  that  an  exact  weak  coupling  calculation  of  the  absorption 
band  shape  would  show  line  mixing  effects  at  least  as  great  as  those  shown 
in  Figs.  1-6  computed  using  the  linear  approximation. 

Several  areas  of  additional  research  are  suggested  by  the  results  of 
this  study.  First,  the  multi-line  Q-branch  spectrum  should  be  rigorously 
diagonalized  in  the  weak-coupling  approximation.  This  is  a  difficult 
numerical  problem  but  the  QR  algorithm^,  or  its  equivalent,  could  be  used 
to  obtain  the  spectrum.  Second,  experimental  absorption  spectra  of  C0o  arc 

t- 

needed  at  pressures  above  10  torr  where,  from  the  results  of  this  work,  line 
mixing  effects  should  be  observable.  It  would,  of  course,  also  be  of  interest 
to  extend  the  measurements  to  pressures  of  several  atmospheres  where  motional 
narrowing  should  be  observed.  Finally,  the  validity  of  the  weak  coupling 
approximation  for  CO2  has  not  been  tested  by  the  present  work.  A  detailed 
calculation  of  the  absorption  band  shape  is  clearly  needed  in  which  elements 
of  the  cross  section  matrix,  given  by  Eq.  (3),  are  computed  directly  from  the 


collision  dynamics. 
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Appendix  A.  Perturbation  Theory  Calculation  of  Eigenvalues  and  Eigenvectors 

of  the  Matrix  (oj  +  i!i) 

0 


Consider  the  eigenvalue  equation: 

(A  +  eB)x  =  Xx  , 


(A.l) 


where  c  is  a  smallness  parameter.  We  expand  the  eigenvalues  and  eigenvectors 
in  a  power  series  in  £  and  collect  terms  proportional  to  each  power.  To 
second  order: 

Ax<°)  =  A(°>x(o)  ,  (A.l) 

Ax(l)  +  Bx(°)  =  X^M1)  +  ,  (A. 3) 

Ax^2^  +  =  A^^x^2)  +  X^^x^^  +  X^2^x^®^  .  (A.l) 


We  assume  in  what  follows  that  the  eigenvalues  A'U) 

k 

the  corresponding  eigenvectors  are  orthogonal 

k 


■  v 


are  non-degenerate  and 
and  normalized  so  that: 


(A 


we  then  expand  x 


O) 


in  terms  of  the  x^°) : 


\u)  ■  I  ■«> 


(A. 6) 


and  it  may  be  shown  that  the  first  order  corrections  to  the  eigenvalue  and 
eigenvector  are: 


(A.  7) 
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a(l>  -  — Si -  ,  t  4  k 

x(»)  -  x<»> 

k  £ 


Similarly,  we  expand  in  terms  of  the  x)0'  : 

k  £ 


x(2)  =  y  (2)  x(o) 

’  1  ik  £ 


(A.  9) 


k  l  2. 

and  it  may  be  shown  that  the  second  order  corrections  to  the  eigenvalue  and 


eigenvector  are: 


x(2)  =  y 

k  £^k 

K.  *■» 


B„  B  , 

(2)  _  r  _ im  mk _ 

5,k  (A(0)_a(°))  (A(o>_xCo)) 

k  £  km 


In  our  case  we  define: 


,(2)  =  r  V&k 
k  ik  V“* 


(A. 10) 


(A. 11) 


\u  ■  <“k +  10k)stk  • 

(A. 12) 

V  -  ‘V  •  :  * k 

(A. 13) 

=0  ,  £  =  k. 

the  eigenvalues  and  eigenvectors : 

(0) 

Xk  =  \  +  iQk  * 

(A. 14) 

=  °* 

(A. 15) 

(A. 16) 
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(A. 17; 


x(1)  =  i 

«k 


Ik 


V 


l  t  k 


(A. 18) 


(?)  - 


UK. 


=  -  l 

m^k 


’’  £m 


(W  -W.  )  (ul  -C;  ) 
k  5,  Iv  m 


1  jt  k 


where  in  computing  A^^  -  A^  ^  =  w^-aj^+i(a^-a^ )  we  have  neglected  the  imaginary 
term  since  the  linewidths  of  adjacent  lines  are  nearly  identical. 
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Appendix  B.  Weak  Coupling  Calculation  of  a  Four  Line  Spectrum 

We  wish  to  compute  the  absorption  band  shape,  F(to),  given  by  Kq .  (11) 


as : 


F(gj)  =  --  Imd*(oil  -  to  -  in)-1*P*d  • 
V  0 


From  Eqs.  (16)  and  (20c),  we  have  for  Q-branch  lines: 


F(m)  =  ±  I  6v[  2  H(u)  , 


where 


H(w)  =  lm(d)l  -  w  -  in)  ^*p 

u 


The  formal  problem  is  then  to  obtain  the  inverse  of  the  matrix: 


(3.1) 


(B .  2) 


(B.3) 


M  =  (wl.  -  u  ~  ill) 
0 


(B.4) 


st 


In  the  weak-coupling  approximation  M  contains  only  diagonal  and  1  off- 
diagonal  terms.  It  thus  has  the  form: 


The  inverse  is 


M 


M 


0 
0 

rl  = 


'll 

mi2 

0 

0 

“21 

m22 

m23 

0 

i 

m32 

m33 

m34 

l 

0 

m43 

m44 

(det  M)-*  M  , 


(  B .  5 ) 


(B.6) 


where 
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det  M  =  b1b2  -  ®11"‘44™23d.32. 


(B.7) 


M  = 


l22b2  ~  m23m32m44 

”mi2b2 

ri12,a2  3m44 

nl2n23n34 

“m21b2 

milb2 

::"1Um44 

mllm23ffi34 

,32m21m44 

“m32mHm44 

m44bl 

"m34bl 

-m43m32m21 

mllm43m32 

~m43bl 

m33brm23m 

In  Eqs.  (B.7),  (B.8)  we  have  used  the  abbreviations: 


(B .  8) 


J 


bl  =  mlim22  "  m!2m21  ’ 


b2  m33m44  ~  m34n43 


( B  .  9 ) 

(B .10) 


We  then  define  the  following  quantities: 


Xk  w_i"k  ’ 


Y2k  V\  "  aZak  ' 


L.,  =  X  a.  +  X,  a.  , 
2k  9.  k  k  I.  ’ 


Z2k  '  Y2k  +  n2knk2  ’ 


(B. 11) 

(B  .12) 

( !i .  L  J) 

(B.14) 


in  terms  of  which  we  may  express  the  real  and  imaginary  parts  of  detM  and  M: 


detM  =  D'  -  iD", 
M  =  ft’  -  ill" , 


(B.15) 

(B.16) 
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In  the  form: 


Z12Z34  +  Y14n23H32  '  L12L34  ’ 


D"  =  212L34  +  Z34L12  +  L14n23I[32  * 


(B  .17) 

(B.18) 


X2Z34  a2L34+II23ri32X4  R12L34 


n21L34 


-n32n21X4 


XlZ34~alL34 


n32L14 


jtI43jTI32Xl 


'Il12n23X4  0 


!!23L14 


X4Z12  a4X12 


H43L12 


■II23X34X1 


j 

1  J.  > 


X3Z12_a3Ll2 


■+•  it  r  x 

23  32  1 


(3.19) 


M’ 


X2L34+'CX2Z34 

+a4n23n32 


'n21Z34 


n32I!21a4 


ri43n32N21 


"  iI12Z34 


alZ34+XlL34 


X32Y14 


n43n32al 


1!12,:23a4 


-n  y 

23  14 


U4Z12+X4L12 


It  Z 
4  3  12 


II12‘!23r34 


-1T  V  a 
23  ‘34  1 


II  7 
34  12 


X3L12+a3/'12 


+  a 


1  23  32 


In  terms  of  these  quantities  the  spectral  function  H(u>)  may  be  written  in 


the  form 


I 


H(ui)  -  (D’2  +  D"2)  1  l  ~  U'Sk1’ 

k  =1 


04. 21) 


where : 


(15.22) 


¥  =  ^ 


(B. 23) 


The  absorption  coefficient  for  the 
by  replacing  Gq(<*»)  in  Lfl*  by  ‘ 


synthetic  4-line  Q-btanch  is  then  obtained 


c0O) 


(B • 24) 


